50.2 SCOPE vs STOP

Fit LME models for Prolactin, HPL and GH2.
Fit GEE models for Insulin and HOMA-IR.

• Outcomes: Prolactin, hPL, GH2, Insulin, HOMA-IR, Glucose (all outcomes log-transformed)
  
• Independent variables: Study, GA at sampling
  
• Model fitted: LME or GEE model
        log(hormones) ~ Study + GA at sampling + Study:GA at sampling
        LME only: Random intercepts for Study ID and Plate
        GEE only: Independence correlation structure

50.2.1 Global tests

Type III ANOVA on main and interaction effects, and LM model summaries.

PRL (ng/ml)

	Chisq	Df	Pr(>Chisq)
(Intercept)	8.706512	1	0.003	**
Study	2.244395	1	0.1	ns
GA_samp	30.790407	1	<0.0001	***
Study:GA_samp	2.221336	1	0.1	ns

hPL (μg/ml)

	Chisq	Df	Pr(>Chisq)
(Intercept)	46.94900	1	<0.0001	***
Study	21.66283	1	<0.0001	***
GA_samp	94.21469	1	<0.0001	***
Study:GA_samp	25.95307	1	<0.0001	***

GH2 (ng/ml)

	Chisq	Df	Pr(>Chisq)
(Intercept)	12.13572	1	0.0005	***
Study	16.86104	1	<0.0001	***
GA_samp	53.71079	1	<0.0001	***
Study:GA_samp	17.66908	1	<0.0001	***

Insulin (uU/ml)

	model term	df1	df2	F.ratio	Chisq	p.value
1	Study	1	Inf	6.440	6.440	0.01	*
3	GA_samp	1	Inf	3.036	3.036	0.08	ns
2	Study:GA_samp	1	Inf	4.677	4.677	0.03	*

HOMA-IR

	model term	df1	df2	F.ratio	Chisq	p.value
1	Study	1	Inf	7.472	7.472	0.006	**
3	GA_samp	1	Inf	3.420	3.420	0.06	ns
2	Study:GA_samp	1	Inf	4.874	4.874	0.03	*

Glucose (mmol/L)

	model term	df1	df2	F.ratio	Chisq	p.value
1	Study	1	Inf	0.232	0.232	0.6	ns
3	GA_samp	1	Inf	0.008	0.008	0.9	ns
2	Study:GA_samp	1	Inf	0.111	0.111	0.7	ns

50.2.2 Model summary

Ratio of Geometric means (95% CI)

PRL (ng/ml)

Var	Effect	P
(Intercept)	1.72 (0.58, 2.87)	0.003
StudySTOP	0.91 (-0.28, 2.10)	0.13
GA_samp	0.21 (0.13, 0.28)	<0.0001
StudySTOP:GA_samp	-0.06 (-0.14, 0.02)	0.14

hPL (μg/ml)

Var	Effect	P
(Intercept)	2.14 (1.52, 2.76)	<0.0001
StudySTOP	-1.51 (-2.15, -0.87)	<0.0001
GA_samp	0.20 (0.16, 0.24)	<0.0001
StudySTOP:GA_samp	0.11 (0.07, 0.15)	<0.0001

GH2 (ng/ml)

Var	Effect	P
(Intercept)	-1.33 (-2.10, -0.56)	0.0005
StudySTOP	-1.64 (-2.45, -0.83)	<0.0001
GA_samp	0.18 (0.13, 0.23)	<0.0001
StudySTOP:GA_samp	0.11 (0.06, 0.17)	<0.0001

Insulin (uU/ml)

Var	Effect	P
(Intercept)	4.31 (2.75, 5.87)	<0.0001
StudySTOP	-2.10 (-3.73, -0.48)	0.01
GA_samp	-0.11 (-0.21, -0.01)	0.04
StudySTOP:GA_samp	0.12 (0.01, 0.23)	0.03

HOMA-IR

Var	Effect	P
(Intercept)	3.01 (1.39, 4.64)	0.0003
StudySTOP	-2.36 (-4.05, -0.67)	0.006
GA_samp	-0.12 (-0.22, -0.01)	0.03
StudySTOP:GA_samp	0.13 (0.01, 0.24)	0.03

Glucose (mmol/L)

Var	Effect	P
(Intercept)	1.67 (1.34, 1.99)	<0.0001
StudySTOP	-0.08 (-0.42, 0.26)	0.63
GA_samp	0.00 (-0.02, 0.02)	0.90
StudySTOP:GA_samp	-0.00 (-0.03, 0.02)	0.74

50.2.3 Diagnostics

Model diagnositcs.

Diagnositcs for PRL (ng/ml)

Diagnositcs for hPL (μg/ml)

Diagnositcs for GH2 (ng/ml)

Diagnositcs for Insulin (uU/ml)

Diagnositcs for HOMA-IR

Diagnositcs for Glucose (mmol/L)

50.2.4 Estimated means

Estimated marginal means (or adjusted means) adjusted for GA at sampling.

PRL (ng/ml)

Study	GA_samp	response	SE	df	lower.CL	upper.CL
SCOPE	10	44.67877	9.765118	775.60673	29.09173	68.61720
STOP	10	60.90187	4.717695	79.27188	52.20009	71.05425
SCOPE	16	155.26985	12.165875	82.42309	132.86150	181.45757
STOP	16	147.84441	14.552006	197.53829	121.76053	179.51605

hPL (μg/ml)

Study	GA_samp	response	SE	df	lower.CL	upper.CL
SCOPE	10	60.06118	6.777997	1013.3294	48.13041	74.94940
STOP	10	39.61791	1.055110	108.7181	37.58090	41.76535
SCOPE	16	194.04839	5.259046	113.7242	183.90473	204.75155
STOP	16	247.31102	10.211412	486.6536	228.03943	268.21126

GH2 (ng/ml)

Study	GA_samp	response	SE	df	lower.CL	upper.CL
SCOPE	10	1.600762	0.2192977	433.4993	1.222895	2.095387
STOP	10	0.948957	0.0264086	216.5576	0.898308	1.002462
SCOPE	16	4.703807	0.1314798	197.9604	4.451542	4.970366
STOP	16	5.462550	0.2814179	922.9232	4.937259	6.043728

Insulin (uU/ml)

Study	GA_samp	response	SE	df	asymp.LCL	asymp.UCL
SCOPE	10	25.31505	7.1937998	Inf	14.504141	44.18404
STOP	10	10.17569	0.4433112	Inf	9.342878	11.08274
SCOPE	16	13.27781	0.7140309	Inf	11.949565	14.75370
STOP	16	10.90728	1.0191172	Inf	9.082069	13.09930

HOMA-IR

Study	GA_samp	response	SE	df	asymp.LCL	asymp.UCL
SCOPE	10	6.371755	1.8874508	Inf	3.565455	11.386841
STOP	10	2.137866	0.0979181	Inf	1.954312	2.338660
SCOPE	16	3.173325	0.1800916	Inf	2.839275	3.546678
STOP	16	2.273611	0.2215450	Inf	1.878337	2.752067

Glucose (mmol/L)

Study	GA_samp	response	SE	df	asymp.LCL	asymp.UCL
SCOPE	10	5.361853	0.3217762	Inf	4.766861	6.031111
STOP	10	4.750210	0.0569715	Inf	4.639850	4.863194
SCOPE	16	5.406086	0.0560985	Inf	5.297246	5.517163
STOP	16	4.682279	0.0697641	Inf	4.547522	4.821030

50.2.5 Comparisons

Hormone	Contrast	GA_samp	Ratio of Geometric means (95% CI)	P-value
PRL (ng/ml)	STOP / SCOPE	10	1.36 (0.90, 2.07)	0.1
PRL (ng/ml)	STOP / SCOPE	16	0.95 (0.81, 1.11)	0.5
hPL (μg/ml)	STOP / SCOPE	10	0.66 (0.53, 0.82)	0.0003
hPL (μg/ml)	STOP / SCOPE	16	1.27 (1.18, 1.38)	<0.0001
GH2 (ng/ml)	STOP / SCOPE	10	0.59 (0.45, 0.78)	0.0002
GH2 (ng/ml)	STOP / SCOPE	16	1.16 (1.04, 1.29)	0.007
Insulin (uU/ml)	STOP / SCOPE	10	0.40 (0.23, 0.71)	0.002
Insulin (uU/ml)	STOP / SCOPE	16	0.82 (0.67, 1.01)	0.07
HOMA-IR	STOP / SCOPE	10	0.34 (0.19, 0.60)	0.0003
HOMA-IR	STOP / SCOPE	16	0.72 (0.57, 0.89)	0.003
Glucose (mmol/L)	STOP / SCOPE	10	0.89 (0.79, 1.00)	0.05
Glucose (mmol/L)	STOP / SCOPE	16	0.87 (0.84, 0.90)	<0.0001

50.2.6 Barplots by Study

Barplots in this section may be misleading due to involvement of Study by GA interaction.

50.2.7 Plots for GA